Purdue 

UNIVERSITY 


Adaptable Single Active Loop Thermal Control 
System (TCS) for Future Space Missions 


Issam Mudawar & Seunghyun Lee 

Purdue University Boiling and Two-Phase Flow Laboratory (PU-BTPFL) 
Mechanical Engineering Building, 585 Purdue Mall 
West Lafayette, IN 47907-2088 

Mohammad Hasan 
NASA Glenn Research Center 
Fluid Physics and Transport Branch 
21000 Brookpark Rd, 

Cleveland, OH 44135 


ASGSR 2015 


Boiling and Two-Phase Flow Laboratory (BTPFL) November 2015 


1 





Research Objectives 

UNIVERSITY 


Condenser/ Compressor 

Heat Out Radiator 



Develop Adaptable Single-Loop Thermal Control System (ASL-TCS) that can tackle 
widely varying heat loads of different space missions and both cold and warm 
environments while improving system performance, reducing weight and volume, and 
ensuring reliable operation irrespective of gravity 
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Thermal Requirement for Different Missions 

UNIVERSITY 


Thermal loads and effective sink temperatures for different missions 
based on Orion Multi-Purpose Crew Vehicle (MPCV) 




TLC, 

TNEOC, 

TMC 

LLO 

LSO 

NEO 

MSO 

LMO 

Thermal 

Load 

1.2 kW 

1 kW 

5 kW 

6.25 kW 

TBD 

6.25 kW 

5 kW 

Effective 

Sink 

Temperature 

-93 to -66°C 
(180 to 207K) 

-198°C 
(75 K) 

-213 to 17°C 
(60 to 290 K) 

-56 to -34°C 
(21 7 to 239K) 

TBD 

-123 to -23°C 
(150 to 250K) 

22°C 
(295 K) 


LEO: 

Low Earth Orbit 

LSO: 

Lunar Surface Operation 

LLO: 

Low Lunar Orbit 

NEO: 

Near-Earth Object 

TLC: 

Trans Lunar Coast 

MSO: 

Mars Surface Operation 

TNEOC: 

TMC: 

Trans Near Earth Objects Coast 
Trans Mars Coast 

LMO: 

Low Mars Orbit 


LLO and LMO associated with heat sink temperatures that exceed lowest minimum 
fluid temperature of 2°C at inlet to cabin HX, therefore requiring heat pump mode 
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T — 0°p- 

^ cabinHX,in ^ ^ ■ 

7 1 — 5°P- 

evap 0 ^ " 


Condenser pinch temp. = 5°C 


Requirements for Heat Pump Mode 


prevent frost formation on evaporator surfaces 

evaporation temperature 

maintain adequate condenser performance 


LLO: T sink = 17°C, T 


cond, out 


- 22°C 


LM O : T sink - 22°C, T condout — 27°C 

Max T lift = 80°C when T evap = 5°C (60°C for R404a with T crit = 72.07°C) 

QcabinHX = 0-75 kW 


Qavionii 


avionicsHX 


= 0-5.5 kW 


Pressure: 


Max pressure = P sat at T c 
Min pressure = P sa , at T cabinHXJn = 2°C 


cond„sat ~ ~^evap 5°C) + T m (— 80°C) 


Radiator Area: 

tit c (t — T ) 

. "'hp ^ p,f\ cond, in cond, sat) 


A =A +A +A T = T +T 

^radiator, s ^ vapor ^ sat t liquid > * cond, sat ± evap ± lift 


vapor 


e a 


T +T \ 

cond ,in cond, sat 

2 ) 


-T 


sin k 


’ A ,a, ~ 




I rri 4 rr,4 

c cond, sat sink 


, A 


m r (T — T ^ 

hp p,f\ cond, sat cond, out ] 


liquid 


£ 0 


\ 


(t +t \ 

cond , sat cond , out 


rri4 
1 sin k 


Coefficient of Performance: COP = 


Q 


W In 

comp / I 


cabinHX Q avionicsHX _ y\ = 

5 I comp 


comp 


comp / i comp 


_ C Pcond - in V dP = m (h -h ) 

K r comp J p v ,u hp \ ‘'cond , in ' 1 avionicsHX , out J 


avionicsHX , out 
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Mass Flow Rate 


Maximum and Minimum Pressure 


0.05 


0.04 


T3 

Eg 

U-2 
o re 
o; o; 

E 5 

3 O 

E lI 

'c </> 
= </> 
S re 


0.03 


ll 0.02 


0.01 


0.00 



Mass Flow Rate: 

- LLO: Lower T condout , lower enthalpy at inlet to 
cabin HX, higher evaporator enthalpy rise, 
lower mass flow rate 


Pressure: 

- Max pressure highest for NH 3 , lowest for 
HFE7000 


- NH 3 requires lowest flow rate 



ASGSR 2015 

Boiling and Two-Phase Flow Laboratory (BTPFL) 

November 2015 

5 




Minimum Required Radiator Area (m 2 ) 
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Low Lunar Orbit (LLO) 


Low Mars Orbit (LMO) 
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Temperature Lift (°C) 


Mass Flow Rate: 

- Area decreases with increasing T lift 

- NH 3 requires smallest area but area differences among fluids are small 










Coefficient of Performance (COP) 
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Coefficient of Performance ( COP) for Heat Pump Mode 


Low Lunar Orbit (LLO) Low Mars Orbit (LMO) 



Temperature Lift (°C) Temperature Lift (°C) 


COP Trends 

- Decreases with increasing T lift because of higher compressor work 

- Higher for LLO because of lower evaporator inlet enthalpy 

- Fairly equal for different fluids 
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Fluid Selection 


1 . 


Chlorofluorocarbons (CFCs) 

• Consist of chlorine, fluorine and carbon 

• Ex: R11 , R12, R113, R114 and R115 

• High ozone depletion potential (ODP) 

• Already banned 

Hydrochlorofluorocarbons (HCFCs) 

• Consist of hydrogen, chlorine, fluorine and carbon 

• Ex: R22, R123, R124 and R142b 

• Small but finite ODP 

• Being gradually phased out 

Hydro fluorocarbons (HFCs) 

• Consist of hydrogen, fluorine and carbon 

• Ex: R404a, R134a, R245fa, HFE7000 

• Zero ODP; low global warming potential (GWP) 

• Currently preferred refrigerants 


Toxicity Rating: 

Class A: 

no evidence of toxicity below 400 ppm 

Class B: 

evidence of toxicity below 400 ppm 

Flammability Rating: 

Class 1 : 

no flame propagation in open air 

Class 2: 

may propagate flame under certain co 
nditions in open air 

Class 3: 

highly flammable 




Type 

Rating 

Critical 

Pressure 

(MPa) 

Critical 

Temperature 

(°C) 

Ammonia 

- 

B2 

11.33 

132.25 

R404a 

HFC 

A1 

3.73 

72.12 

R134a 

HFC 

A1 

4.06 

101.06 

R245fa 

HFC 

B1 

3.65 

154.01 

HFE7000 

HFC 

A1 

2.48 

165.00 

R123 

HCFC 

B1 

3.66 

183.68 
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Comparison of Fluid Performances 

UNIVERSITY 



NH 3 : 


Requires smallest flow rate but highest condenser pressure; exhibits toxicity and 
flammability 


R404a: An HFC that requires high condenser pressure; exhibits no toxicity or 

flammability, zero ODP and low GWP 

R134a : An HFC that provides compromise between reducing flow rate and reducing 

pressure; exhibits no toxicity or flammability, zero ODP and low GWP 

R245fa: An HFC that provides compromise between reducing flow rate and reducing 

pressure; exhibits no flammability, zero ODP and low GWP, but some toxicity 

HFE7000: An HFC that requires low condenser pressure but high flow rate; exhibits no 

toxicity or flammability, zero ODP and low GWP; used as replacement to R123 


R123: 


An HCFC that requires low condenser pressure but high flow rate; exhibits no 
toxicity or flammability but some ODP 


No significant differences among fluids in terms of radiator area or COP 


Preferred Coolant: R134a (1,1,1,2-Tetrafluoroethane, CH 2 FCF 3 ) 
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50 100 150 200 250 300 


h [kJ/kg] 

- LLSL-HX used high temperature of subcooled liquid exiting condenser to further superheat 
vapor before entering compressor 

- Compared to basic heat pump: 8 to 8’ and 4 to 4’ 
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T-S Diagram for Heat Pump with LLSL-HX for Low Mars Orbit (LMO) 

using R134a 
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o 50 

o 


T evap = 5°C, T lift = 50°C 
n C om P = 0-85, e reg =0.9 



0.50 

s [kJ/kg.K] 


COP' = 


( QcabinHX QavionicsHX ) ^Q ( 


cooling 


W + Z1W 

comp comp 



cooling 


QcabinHX Qavii 

AW^„ 


avionicsHX 


COP 
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Heat Pump with LLSL-HX versus Basic Heat Pump Mode for R134a 

UNIVERSITY 
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3.2x10-2 


3.0x10-2 


LMO (Heat Pump) 
LLO (Heat Pump) 


2.8x10-2 



LMO (Heat Pump + LLSL-HX) 
LLO (Heat Pump + LLSL-HX) 



Evaporation Temperature (°C) 



Temperature Lift (°C) 



Temperature Lift (°C) 


Low Lunar Orbit (LLO) 
Low Mars Orbit (LMO) 


Mass Flow Rate: appreciable reduction with 
LLSL-HX 

COP: some degradation with LLSL-HX 

Radiator Area: some improvement with 
LLSL-HX 
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Heat Pump Test Loop 
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Avionics HX: Cross Section and Temperature Measurements 



i 


Type-E, Miniature 
Thermocouple 



©©©©©©©©©© 


Thermocouple 
at inlet 

Thermocouples inside copper block 

Thermocouple 
at outlet 
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Heat Transfer Characteristics 

and 

Dryout Effects 
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Purdue Avionics H/X: Two-Phase Flow Patterns and Heat Transfer Mechanisms 

UNIVERSITY 
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*e 

■ Intermittent Dryout: mostly in liquid film surrounding elongated slug flow bubble or due to vapor blanket 
formation within the liquid slug 

■ Incipient Dryout: dry wall patches in annular film 

Dryout Completion: annular flow replaced by mist flow, with cooling provided only by droplet deposition 
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Purdue Avionics H/X: Two-Phase Flow Patterns and Heat Transfer Mechanisms 

UNIVERSITY 


(i) x e < 0.36 



Liquid film Intermittent dryout due Bubble formed Intermittent 
to vapor blanket by liquid dryout in 

formation in liquid slug bridging liquid film 


Annular flow 



Wavy annular liquid Liquid droplets 

film shear driven by entrained in vapor 

vapor core core 




Annular flow 

► 


(iii) 0.50 <x e < 0.74 



Incipient dryout caused Dry wall patches 

by formation of dry cooled only by 

patches in liquid film droplet deposition 




Purdue Assessment of Heat Transfer Coefficient Correlations 

UNIVERSITY 


Macro-channel: 

Chen (1966) 

Liu & Winterton (1961) 

Mini/Micro-channel: 

Lazarek & Black (1982) 
Tran et al. (1996) 
Kandlikar (2004) 

Lee & Mudawar (2005) 
Bertsch et al. (2009) 
Kim & Mudawar (2014) 
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Comparison of Present Average Two-Phase Heat Transfer 
Coefficient with Predictions of Universal Methodology 
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Universal Predictive Methodology for Saturated Flow Boiling 
Heat Transfer in Small Diameter Tubes 



- Working fluid: 

FC72, R11, R113, R123, R1234yf, R1234ze, R134a, 
R152a, R22, R236fa, R245fa, R32, R404A, R407C, 
R410A, R417A, C0 2 , water 

- Hydraulic diameter: 

0.19 <D h < 6.5 mm 

- Mass velocity: 

19 < G< 1608 kg/m 2 s 

- Liquid-only Reynolds number: 

57 < Re f0 < 49,820 

- Flow quality: 

0 < x< 1 

- Reduced pressure: 

0.005 < P R < 0.69 

Kim & Mudawar, Int. J. Heat Mass Transfer 64 
(2013) 1239-1256 


For nucleate boiling dominant regime : 


h nb = 2345 Bo - JL 


P R 038 (l - x e )"°' 51 0.023 Ref Prf 


For convective boiling dominant regime : 

( p r 8 ( i \ o94 (p 

h,= 5.2 Bo^] We,-~°' 54 + 3.5 — 


x 0.023 Re 0 / Pr° A — 

i / n 


Q 

where Bo = - UL - 
Gh, 


, P R = 


Re f = 


G(l~x e )P h 

ix f 


We /o = 


( u rvi-r \ o9 (p r 
^ p ^ 

\ x„ p e 


n 

Qh. effective heat flux averaged over heated perimeter of channel 
P H : heated perimeter of channel 
P F \ wetted perimeter of channel 
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Universal Predictive Methodology for Saturated Flow Boiling 
Heat Transfer in Small Diameter Tubes 


Nucleate Boiling Dominant Data 



0.5 1 10 100 

h tB (exp) [kW/m 2 K] 


Dryout Dryout 
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I 1 1 



Nucleate Boiling Dominant Heat Transfer 


Convective Boiling Dominant Data 



Dryout Dryout 

= 0 incipience completion 
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G crew = 264.63 kg/m 2 s 
q”c r ew = 16146.9 W/m 2 
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Small amplitude, high frequency pressure drop oscillations 
(4,000 - 1 0,000 Pa over period of 1 .33 s ) 
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Averaged Power Spectrum [Pa 2 ] 
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pressure drop response in frequency domain 
(power spectrum peak of 0.74 Hz (~ 1 .33 sec) 
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